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^ ", Abstract 

o 

^ ■ We investigate coupled gravity and a bulk scalar field on a slice of AdSs bulk 

. with special BPS branes at the two ends. With the special scalar potentials on the 

Q \ branes the scalar field does not stabilize the size of the orbifold. With a careful 

l/^ ' treatment of the general coordinate invariance the complete tensor-scalar spectrum 

\^ . is presented. There are two massless zero modes in the scalar sector, the radion 

\ and a dilatonic zero mode. The scalar KK modes acquire masses at the order of the 

' warped mass scale. The four dimensional effective action is of tensor-scalar type. 

> ' 
(N , 

^ 1 Introduction 



J3 



■ Recently there has been a growing interest in extra dimensional models where the 

, standard model fields localized on a brane. Randall and Sundrum proposed a sce- 

^? \ nario (RSI) with one extra dimension, two branes of opposite tension at the orb- 

' ifold points of an AdS bulk [1]. This setup was motivated by the recent progress 

^ • in string/M-theory and can be realized in string theory [2], [3]. Poincare invariant 

D , solution requires a finetuning of the bulk cosmological constant A = — and the 

'■^ brane tensions Xpianck = —^TeV = The exponentially warped metric can gener- 

ate the hierarchy between the Planck and the electroweak scale with moderate value 
^ ■ parameters krc ~ 0(40). However, there is an exact solution for any Vc distance of 

the two branes. A change in Tc is described by a four dimensional massless scalar 
field, the radion [4]. It couples to matter like a Brans-Dicke scalar and must have 
a large Brans-Dicke parameter or a mass to recover standard 4-dimensional gravity 
[5], [6]. 

Goldberger and Wise proposed a mechanism with a massive bulk scalar field 
to stabilize the radion [7] . They added special steep Xi{^) brane potentials. The 
competition between the kinetic and potential terms generates a minimum in the 
effective potential and determines Vc- Approaches based on naive ansatz without 
including the backreaction found that the radion acquires an O(TeV) mass and the 
couplings to standard model fields are ~ 1/TeV [5]. Effective field theory calculation 
gave the same radion mass in [8]. Tanaka, Montes [9] and later independently 
Csaki et al. included the backreaction and derived a single equations for the scalar 
degrees of freedom based on Lorentz covariance of the linear perturbations [10]. The 
equation together with the boundary conditions at the branes defines a hermitic 
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operator only in special cases, like the stiff potential limit, when — ^ oo. In this 
limit the physical results agree with the naive calculations. 

Another mathematically well defined case is when the brane potential is equiva- 
lent to the superpotential that generates the bulk scalar potential (13), we call these 
BPS brancs. This type of models with exponential potentials can possibly arise from 
the bosonic sector of dimensionally reduced supergravity theory [11]. In this letter 
we give the complete spectrum of linear perturbations around an exact background 
solution of two flat BPS branes in 5 dimensions. This scenario was analyzed in 
[12] using two different Gaussian normal coordinate patches with respect to the two 
branes following the method of [4]. They discussed only one scalar perturbation 
the massless radion. Solving the equations with two flat branes we find still one 
more massless dilatonic scalar and the KK tower for the scalar field. This model 
does not solve the moduli problem but the complete pcrturbative spectrum can be 
understood, the limit to the Randall-Sundrum case is smooth and we hope that it 
gives help to solve the general case with detuned brane potentials. The low energy 
effective theory is of scalar-tensor nature and generally it cannot be accepted by 
phenomenology and cosmology. 

The paper is organized as follows. In section 2. we present the setup and discuss 
the background solution using the technique of [13]. In section 3. we consider the 
simple limit of the Randall Sundrum scenario with an additional bulk scalar field. In 
section 4. we fix the gauge, solve the linearized equations and present the complete 
spectrum and we conclude in the last section. 



2 The model 

In a 5 dimensional ( = {x'^jy = x^) ) M4 x S1/Z2 manifold, there are scalar 
field, gravity and two flat branes at the orbifold fixed points y = 0, Vc- Latin indices 
run from 1..5 and greek indices 1..4 and prime denotes the derivation with respect 
to the proper coordinate y. 4-dimensional greek indices raised and lowered by r/^,/. 
The action we consider is 




Qi is the induced metric on the brane, = 1/(2M^) is the 5D cosmological constant, 
the integration understood on the orbifold. 

We want to investigate perturbations around a background solution. The ansatz 
for solving the Einstein scalar coupled system that respects Lorentz invariance is 

ds'^ = a\y)ri^^dx''dx'' -dy^, (2) 
$(x,y) = $o(y)- (3) 

The warp factor a(y) and $o(y) are continuous even functions. The first derivatives 
may have jumps, the second derivatives delta singularities at the branes. We use 
the following form of the Einstein equations 

Rab = l^'^fab = (Tab - Igabg'^T^d] ■ (4) 
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The gravitational and scalar equations are the following. 



fa'V fa'\' , «^v-. 



3 
a 



(5) 
(6) 
(7) 



The boundary conditions are 



-I Vi 



A- 



2/« 



(8) 



With Z2 symmetry and working in the (0, Vc) interval the jump is twice the value of 
the derivative calculated from the inside of the interval and there is sign difference 

at the two branes. 

An analytic solution of these equation is generated by the technique proposed 
by DeWolfe et al.[13]. With one superpotential a flat solution is generated 

via 



y($o) = 



2 5$ ' 

-^VF($o). 

D 

1 fdW{^Q) 
8 V 



(9) 
(10) 

(11) 



These solve the equations of motions if the boundary conditions for 1/F(<&) are ful- 



filled (+ for 2/ = and - for y 



-- ±W{^^) 



(12) 



Exactly one fine tuning is needed to find a flat background solution. There is 
no need for supersymmetry to find a solution though the bulk plus brane action 
completely vanishes substituting the solution. This technique reduces the coupled 
system of differential equation to decoupled ordinary differential equations and can 
be solved (generally numerically) for arbitrary V{^) and Ai($) tuning one additive 
constant in A2($). Explicit solution can be constructed defining first W{^), like 
W{^) = c - or VF($) = e"*. 

We will investigate perturbations in models where the fine tuning of the scalar 
brane- and superpotential (12) is valid at all order for arbitrary 



Ai($) = ±W{^) 



(13) 



Wc call this type of brancs BPS branes. These branes allow solutions in the bosonic 
sector respecting 5 dimensional bulk+branc supersymmetry, the BPS condition pre- 
serves half of the supercharges corresponding to unbroken N=l supersymmetry in 4 
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dimensions [15]. This type of models with exponential potentials can possibly arise 
from the bosonic sector of dimensionally reduced supergravity theory [11]. The 
localization of scalar modes were discussed with one brane in [16] based on 0(3) 
symmetry in the generalized longitudinal gauge of [17]. 

There is always a flat background solution for any value of Vc- The change of 
Tc is described by a massless modulus field, the radion in the low energy effective 
theory. There is an additional dilatational invariance of the action with BPS branes 
leading to a second zero mode. The value of the scalar field on the Planck brane 
is an integration constant, a free parameter of the action. Choosing a definite 
background solution this symmetry is broken spontaneously and the corresponding 
massless Goldstone boson is expected in the scalar spectrum. In the limit of the 
Randall-Sundrum scenario this invariance is obvious and is presented in the next 
section. 

To obtain an analytical solution for the KK modes the following particular po- 
tential will be used 

(ii. 

W{^) = ±Ai($) = ^e"*. (14) 
The solution contains two integration constants cq, Cp, usually both chosen to 1 

a{y) = co{Cp ^ky)^, (15) 

g-a4.o(2/) = ^Cp-^ky). (16) 

In the limit a ^ we get back the AdS e"'^^ profile of the Randall-Sundrum 
scenario and the scalar field decouples from gravity. We deal with solution with 
< < to avoid naked singularities in the bulk. This solution generates a 

2 

hierarchy of a(rc) = (1 — ^^krc)sa^ . The dilatational invariance corresponds to 
the free choice of Cp = e""*"^*^^ the value of the scalar field on the Planck brane. 
It's change can be compensated by a change in k,co to have the same warp factor 
a{y) and superpotential. For particular solutions we use cq = Cp = 1. In the next 
section we discuss the special limit of the Randall-Sundrum scenario. 



3 RSI with bulk scalar field 

In the RSI scenario specially tuned bulk and brane cosmological constants provide 
a warped solution. In a more realistic model a dynamical mechanism is needed to 

produce the brane and bulk constants. A simple proposal is to add a bulk scalar field 
which has constant bulk and brane potentials. This is the limit of the previously 
presented model with 

W{^) = A($) = ^ = constant, (17) 
y($) = A = - — . (18) 
The classical background solution 

a(y) = e"*''^!, $o(y) = = constant (19) 
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has two invariances. One is related to the free choice of the distance between the 
two branes, the other is the freedom that we can shift the constant scalar field ^rs 
to + ci, it is also an exact solution with the same warp factor. 

In the linearized theory the former provides the massless radion [4] the latter 
zero mode can be identified as follows. The linearized gravitational (hab) and scalar 
{(p{x,y)) perturbations around (19) decouple, it was also realized in [16], [20]. The 
gravitational perturbations are as in RSI, a massless radion with the wavefunction 

F{x, y) = ~^f{x) = k^fix), □/(x) = 0, (20) 

a massless spin-2 graviton localized on the Planck brane and its KK tower. 
The scalar field perturbations are determined from 

D5ct>ix,y)=0, (21) 

with the boundary condition (l)'{x, yi) = 0. This is the free Klein-Gordon equations 
on an AdSs background, with □ = 04 = d^^d^ 

U^ = -dl-A'^d, + ^U. (22) 

The scalar field spectrum contains a dilaton zero mode 

4>{x,y) = h{x), □/i(x)=0, (23) 

reflecting the invariance of the action against the shift of the constant value of 
the scalar fleld. The background solution brakes this symmetry spontaneously and 
this zero mode is the corresponding Goldstone boson. This zero mode mixes with 
gravitational perturbations if the linearized scalar equations are coupled to gravity. 

We solve the equation for the KK modes in Schrodinger basis: change to conform 
flat background by dy = a{z)dz with a{z) = 1/(1 + k\z\) and rescale the fleld by 
(f>{x,z) = e'P^$Qa~20(2;) also changing to mass eigenvectors with D^j = — mf^j. 




The boundary condition is (a~2^(^))' = 0. The potential term multiplying ^ can 
be written as 

The massive solutions are determined by Bessel functions 

4> = ^/l + Wl {aJ2 (^(1 + k\z\)) + BY2 (^(1 + k\z\))) . (26) 

The ratio A/B determined from one boundary condition the mass eigenvalues are 
determined at the other brane and are of the order of the TeV scale. 

The scalar spectrum then contains two zero modes and a KK tower. Turning 
on smoothly the bulk and brane scalar fleld potentials continuous deformation of 
the spectrum is expected. Quantum corrections are also expected to develop a 
nontrivial bulk and brane potential for the scalar field. The linearized equations for 
the action (1) is presented in the next section. 
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4 Linear perturbations 

In this section we give the equation of motion for general Hnear perturbations keep- 
ing the branes at the fixed orbifold points and discuss the role of gauge freedom. A 
general perturbation is 

SQimu = a^hf,„{x,y), Sg^^ = h^v {x,y) 
S955 = = -2G{x,y), 

d^ = c^ix,y). (27) 

Infinitesimal general coordinate transformations 

+ + (28) 

must fulfill the orbifold symmetry conditions, continuous and even, ^5 odd and 
vanishes at the branes. 

C,i{x,y + yi) = ^f,{x,y - yi) (29) 
C5{x,y + yi) = -C5{x,y - yi) (30) 

^5 gauge transformations moving the orbifold fixed points are not allowed in this 
approach. The infinitesimal change of the background solution is understood as a 
change in the perturbations. 

= i^n,u + Cm) - 2^?5??M'^ 

The offdiagonal h^^ components can be transformed out with 

^f^ = - J ^y'Jiy)^^!'' = ^- 

The branes are not moved in the 5*^ dimension. 

Without loss of generality the Einstein equation can be used without offdiagonal 
linear perturbation. The metric is 

ds^ = a^{y){T]^, + h^,{x,y))dx^'dx'' - (1 + 2G{x,y))dz^ . (32) 

(31) was our first gauge choice and there is a restricted gauge degree of freedom to 
keep = . ^5{x,y) is constrained by the orbifold conditions (30) and 

ry I 

^iJi{y,x) = - dyi M A^^yi) + ^i^i^)- (33) 
Jo o (yij 

In what follows we work on the interval (0,rc) and get the fields on the orbifold 
using the Z2 symmetry. 

The linearized equations are the following, '55', '5//' and '//z^' respectively. 



«2 
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3^G^ + I Ka - h,y = ^^'o{y)4>A^,y) (35) 



a'aG' + 2 ( 4a'2 + a2 ( - V ) G 



1 /Alt \' -^^7/2 

+\ [^hl,.)X - - ^.^'^) = (36) 

Here h = h'^, U = d'^d^. The scalar equation is 
1 , ,,, a', dV , 

(5*' - «') *i - - E (0* + °W - y-'- p^' 

The singular terms at the branes define boundary conditions. The complete 5/x 
equation is valid also at the branes. The 55 boundary condition is the trace of 
liu boundary condition calculated from (36). It is equivalent to the Israel jump 
equation calculated with the extrinsic curvature . From /iv (36) we get 

a' 

- 2— G(x, y)r]^u + —<^'o<f>Vtiu = at y = y^. (38) 

Qi O 

Already the leading order boundary conditions and Z2 symmetry were used. The 
scalar equation gives 

±20'(x,yO = ^Hx,yi) + ^Gix,yi). (39) 

The equations are not idependent. The Bianchi identities ensure that from the 
15 Einstein equations are only 10 independent and the bulk scalar equation can 
be derived from the Einstein equations with differentiation and taking linear com- 
binations. The scalar boundary condition presents a meaningful constraint. The 
equation are solved for special BPS branes fulfilling (13). 

4.1 Lorentz decomposition of h^iy 

To solve the equations we decompose h^i, in a Lorentz covariant way. It contains 
a TT (transverse traceless) symmetric tensor (5 components), one divergencefree 
vector (3) and two scalars (1+1) giving 10 components. 

h^iu = hfj,y{x,y) + A^,Ax,y) + Ay^x^y) + 5^5j.6(x,y) - F{x,y)r]^,y (40) 

Here 

K,- = ^. H = o, (41) 

= 0- (42) 
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The components can be uniquely defined if the perturbations are assumed to be 
Lorentz covariant as the background solution is Lorentz invariant and fulfill 

□A^T^O, 06/0, (43) 

to have nonvanishing trace or divergence in h^i,. They are defined by 



n'b = 3 ( - DT^a/? ) /i"^, (44) 

F = ^(-h + ab), (45) 

= (^h\^^ + F^-ab,^). (46) 

This sequential definition of components ensures that F, Ob describe the complete 
trace and F, b, OAfj, describe the whole divergence of h^i,.^ The TT part is defined 
as the vector and scalar parts subtracted from hn^. h^^y can still contain scalar and 
vector TT combinations but those will be proved to be gauge degrees of freedom as 
in four dimensional gravity. 

Gauge transformations of the components can be defined as 

H = -foe5(x,y), (47) 
8G = -6(^,y)', (48) 

SF = 2^^,ix,y). (49) 

The last two transformation can be used to fix ^5(0;, y) without moving the branes 
and setting the F{x, y) — G{x, y) difference only x dependent. Make a gauge trans- 
formation with 

a'^5(x,y) = - r a\yi){Fix,yi)-G{x,yi))dyi+x{y) r a\F-G)dyi, 
Jo Jo 



x{y) = 



fo" a'^{yi)dyi' 



This ^5 fulfills the orbifold condition and it is completely fixed. The gauge transfor- 
mation was motivated by the work of Kubyshin et al. who have found the radion 
wave function with fixed flat branes [14] . In the new gauge 

F"(x, y) - G^{x, y) = /a(x) = ^ a\F - G)dy^ (50) 

Jo o- Ky\m\ Jo 

We work in this gauge and drop the n index in what follows. The remaining gauge 
degree of freedom is 

= e>^{x), (51) 

which can be only used to set the two helicity state of the massless four dimensional 
graviton. The gauge transformation of the other components is 

nSb = 2e^^ ,ifne^^/0 (52) 
6A^ = , if De^ / 0, e^^^ = (53) 



■^A subtlety of the decomposition that h^^ only defines 0^6 while in the definition of F, purely Uh 
appears. 
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Untill this point the equations of motion were not used and the gauge could 
be partially fixed. In what follows we use the Einstein equations and show that 
Afj_{x, y), b{x, y) can be gauged away and h/^i, contains only the graviton and its KK 
tower. The single F{x,y) field will describe the massless radion and the scalar field 
together. 

4.2 Solving the Einstein equations 

The equations can be written up with the general decomposition of h/j^i^. The 5/x 
equation (35) gives 

Ida; = K'^'o{y)^A^, y) - ^F', - ^. (54) 
Different Lorentz spin perturbations decouple and we get 

A'^{x,y)=Q, (55) 
as it was defined as ^ 0. The solution is 

A^^ix, y) = Uix), f% ^ 0, Uf^{x) + 0. (56) 
It can be gauged to nothing with 

= -f^{x). (57) 
The RHS of (54) gives setting one y-dependent integration constant zero 

^a2$'o<^ = ia^G)'. (58) 

In the fiv Einstein equation (36) as the result of assumed Lorentz covariance the 
d^di, derivatives terms decouple from the TT and rjfj^i, part and give the equation 

l(F-G),^. + ^(a%,)' = 0. (59) 
The general solution sending two integration constant in the infinite x to zero is 

b{x,y) = bi{x) + b2{x)J^ + 2/a(x) J ^Ja\if □/a(x) ^ 0. (60) 
The trace of the /iv boundary condition (38) reads 

Db' -4(^F' + 2^G - ^aH'oci^ = 0. (61) 

The (58) solution of 5/x is also valid at the branes giving with Db ^ 

b'{x,yi) = 0. (62) 

At y = it sets 62(2;) = and at rc /a (a;) = 0. The remaining solution in b{x,y) 
can be gauged to zero by 

e\{x) = -lb,{x). (63) 
9 



If ^fA{x) = then using the TT part of the '/Ltz/' equation, one gets f/\{x)^^jj = 
and /a (a;) = 0. 

We proved using the Lorentz covariance of the Unear perturbations that the 
following gauge can always be chosen. There is a TT tensor and a single F{x,y) 
scalar perturbation {F = G). 

h^^ - F{x, y)rji^u (64) 
-2F{x,y) (65) 

F'{x,y) + 2-F{x,y) (66) 
a 

This gauge is called Newton gauge by Tanaka and Montes in [9] and was used in 
[10]. To our knowledge this gauge choice is correct in the presence of two BPS 
branes. The remaining gauge degree of freedom, □e^(a;) = can transform out the 
3 non-physical component of the massless graviton. 

The scalar and the spin-2 TT perturbation decouple. The TT equation is 

^{a%J)' -^ahj:j = 0, (67) 

together with the boundary condition 

h'i,,ix,yi) = 0. (68) 

The boundary condition only allow the massless solution of 

V = h%ix) (69) 

the massless graviton with 2 helicity states. It is localized at the Planck brane with 
larger warp factor. There is also a massive graviton KK tower, with 5 degrees of 
freedom. The KK modes are not localized. 

Upon using the Newton gauge relation (66) the Einstein boundary conditions 
automatically fulfilled. In the bulk the scalar part of the Einstein and the Klein- 
Gordon equations are equvalent. 

This is accompanied by the scalar boundary condition (39). Eliminating <p{x,y) it 
is 

This boundary condition contains a second derivative, it makes generally the opera- 
tor in the radion differential equation non-hcrmitic and causes the hermiticity prob- 
lem observed in [10]. The different mass eigenmodes are not orthogonal to each other 
and does not serve as a basis for expanding a general solution. The bulk equation 
can be used to eliminate the second derivative but then the d'Alcmbcrtian appears 
with different eigenvalues for different eigenfunctions. This simplified boundary 
condition is 



hfxu — 

^55 = 
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It is generally still not hermitic, though with BPS branes the first term vanishes as 

nF{x,yi) = (73) 

boundary condition defines a hermitic operator for massive perturbations. Practi- 
cally there is no boundary condition for the massless solutions. 

4.3 Massless scalar perturbations 

There are two independent massless modes as the boundary condition (73) is auto- 
matically fulfilled. The radion differential equation can be written as 



1 



ii) ar 

The wavefunction of the massless radion is 




(74) 



FR{x,y) = FR{y)r{x) = e{y) irr^r{x) Ur{x) = 0, (75) 

where e(y) is the sign function and it makes the radion wavefunction even under 
Z2. The radion appears in the scalar field perturbation as (j)ji{x,y) = ~2^f{x). 
In the limit of no scalar field (75) gives back the result of Charmousis et al. [4] 
Fn = —ke'^^^r{x). The radion describes the relative free motion of the branes. 
With a gauge transformation of 

^5ix,y) = -^r{x), (76) 

the radion solution is completely transformed to a scalar TT mode and the relative 
motion of the branes. Keeping h^^ = with (33) after this gauge transformation 
the bulk metric is 

f-y 1 

h^iu = Kv + r,^„{x) j Fs{x,y)r]^^ (77) 

Jo o, 

h55 = -2F,{x,y) (78) 

^<^'o^s = F',ix,y) + 2-Fs{x,y). (79) 

Fs,(j)s here only describes the scalar field perturbations. With BPS branes this 
gauge transformation is well defined and the relative motion of the branes is 

This single radion solution was analyzed in [12] without the dynamics of the scalar 
field. 

The scalar field has a zero mode and a complete KK tower. The even massless 
scalar solution of (74) with flat branes is 



a a^{y) Jo 



FD{x,y) = FD{y)fiix) = ^hix) / dy^a'^ , □/i(x) = 0. (81) 

(t) 
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This is the dilatonic zero mode of the scalar field, it is the Goldstone boson of the 
spontaneously broken dilatational invariance of the action. It appeared as adjusting 
a meaningless constant in the limit of the Randall- Sundrum scenario with decoupled 
scalar field. Brax et al. did not discuss this second zero mode in [12]. They were 
aware of two scalar zero modes [18] and recently studied the cosmological evolution 
of them [19]. 

The extra dilatonic zero mode given in (81) vanishes on the Planck brane and its 
wavefunction is the radion wavefunction nuiltiplicd with an extra y-function. The 
two solutions are not orthogonal to each other, they can be orthogonalizcd with 
Gram-Schmidt method, once a scalar product is defined on the interval (0, Tc). Two 
solutions can be distinguished also on physical grounds. The dilatonic zero mode 
(81) also describes a relative motion of the branes through changing the proper 
distance J^" \/\g55\dy = Vc + j^" F{x,y)dy. A new dilaton wavefunction can be 
formed with linear combination that keeps the proper distance Vc unchanged. 

Fn(v) 

F2{y) = Fn{y) - fr^Mv) (82) 

Jo ^Riy) 

Now the two masslcss modes have different physical origin as one moves the branes 
the other not. The radion orthogonal to F2{y) with respect to a specific scalar 
product can be defined as Fr(j/) = Fji{y) - F2iy)iFji, F2)/ {F2, F2). 

The dilatonic zero mode remarkably simplify with exponential superpotential. 
Performing the integral in the solution (81) we get 

FD{x,y) = -h{x)^^ (l^^ + 1) ^ fi{x){FR{y) - F«(0)). (83) 

Subtracting the radion part we get Fj:){x,y) = fi{x)/{l + |^) which implies via 
(66) the dilaton scalar field solution 

cl>i,{x,y) = ~-^Mx). (84) 

In the limit a ^ it is the dilaton of (23). This solution can be found directly as 
the coefficient of the term F in the radion equation (70) vanishes for W ~ e"*, 

f^)'] = 5!,,- nil)' = 0, here W = 55^. (85) 

Thus the equation (70) has a y independent Fuix^y) ~ fi{x) dilaton solution of 
(84). 

The dilatonic zero mode, called 002 was found by Bozza et al. [16] with one brane 
in the case of dilatonic superpotential. This was the only propagating massless mode 
as the radion uj^ is not normalizable and not physical with one brane. UJ2 was iden- 
tified as the dilaton in the limit of the RSI scenario when the scalar field decouples 
from gravity. The Einstein equations gave constraints between the massive modes 
and there were two independent scalar KK tower. One is the scalar field KK tower. 
The other completes the two component transverse traceless tensor hij with two 
additional 3-vector mode to a 5 componenent massive 4-dimensional KK graviton. 
Bozza et al. investigated the scalar perturbation of a Lorentz invariant setup based 
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on 3 dimensional spatial covariancc. At the end 4 dimensional covariant equations 
described the perturbations. It would be enlightening to perform their analysis in 4 
dimensional language and understand the physics of perturbations. Indeed follow- 
ing their method with two branes we also got two massless scalar modes. This work 
is in progress. 



4.4 Massive KK tower 

There is a KK tower of the massive scalar field perturbations with the well defined 
hermitic boundary conditions F(x, yi) =0. In a scenario with small backreaction 
like in [10], the KK spectrum can be solved with good approximation using the 
exponential warp factor a{y) = e~^y . The mass eigenvalues determined by the 
zeros of the Bcsscl functions and are at the order of the weak scale {ke^^^" ). 

The radion equation (70) for the case of exponential superpotential (14) can 
be analytically solved by Bessel functions. As a standard technique we change 
to conform flat background by dy = a{z)dz and rescale the field F by F{x,z) = 
e^P^^Qa^2 F(z) also changing to mass eigenvectors with DFj = —mfFi. The radion 
equation in this Schrodinger basis is 

[4\aJ 2\a) a ^'^ \%) ^ ' 

This equation for each eigenmode has a mass parameter and two integration con- 
stants. One integration constant is fixed by the normalization, the other from the 
boundary condition at the Planck brane {y = 0) and the mass is determined from 
the TeV brane (rc) boundary condition. With the exponential superpotential the 
equation is 



where 



fe(l-^) k 
1 3 



+ z = l + z, (88) 



2 2 _ l' 



(89) 



The range ofi/isl<i/<oo,i/ = l corresponding to the RS scenario. The solution 
fulfilling the y = z = boundary condition is given by two type of Bessel functions 



F=Vc(yu ( Ju (mC) - J. (mo) 



(90) 



The mass eigenvalues arc determined at y = rc from the equation 



with 



-^•"{^ I 7 — r " V^iCc) = 0, (91) 



1,, 3a^ , ,jv2 a(rr) 

{I- krc)^ = (92) 
k 1^ k 
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There is always an mj = eigenvalue, a massless solution. When the warp factor 
generates a large hierarchy between the two branes a(rc) <C o(0) = 1 then the lowest 
mass eigenvalues are practically determined from the second term, the zeros of the 
Bessel function {mi(c)- The lowest KK masses are few time the warped mass 
scale ka^^{rc), the TeV scale. 

F uniquely defines (p via (66) and 

In the limit a = 0, 1/ = 1 we get back the masive spectrum given in (26). 

5 Conclusion 

In this paper we have analyzed the coupled scalar gravitational system in 5 di- 
mensions. Around an exact background solution we carefully fixed the gauge and 
solved the 5 dimensional linear equations of motion on the orbifold. We have found 
that with the special BPS brane potentials there are two scalar zero modes and 
one KK tower. The two scalar zero modes are the radion and the dilaton which 
can be identified in the smooth limit of the RSI scenario with a free bulk scalar 
field. Both modes couple to the matter on the TeV brane via the induced metric. 
The two massless scalar modes generally lead to unacceptable phenomenology, the 
radius is not stabilized. Cosmological evolution of two scalar moduli in this model 
were discussed recently [19]. However the BPS brane scenario may help us to find 
the correct linear spectrum in the Goldberger-Wise stabilization mechanism. The 
gauge fixing discussed in this paper and used in [10] [9] leads to non-hermitic dif- 
ferential operator for scalar perturbations for general brane potentials. The GW 
bulk scalar with detuned potential gives mass not just to the radion but also to the 
dilaton mode. It would be very interesting to follow the fate of these two modes and 
relating them to the first scalar KK mode. Effective theory and AdS/CFT based 
calculation so far could only capture the radion mode [8], [21], [22]. The possible 
resolution of the hermiticity problem is that a mode is droped during solving the 
Einstein equations which has the the detuning of the brane potentials as a source. 
This work is in progress. 
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